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used. To obtain the loads, a three step process was utilized. First, from aircraft 
differential pressure transducers and accelerometers, the spectral content of the response 
and excitation was determined Then, using a NASTRAN model modified to replicate the 
flight test aircraft, frequency response functions were determined between selected points 
on the aircraft's tail and the accelerometers. Finally, using this information, a solution was 
obtained for the vibratory airloads. Having provided information on the nature of the 
driving forces, structural modifications can be made that move the natural frequencies 
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I. INTRODUCTION 


A. GENERAL 


The RAH-66 Comanche is the United States Army's newest light/attack armed 
reconnaissance helicopter. It is designed to be a combat multiplier on the joint and 
combined arms battlefields of the 21“ century. The aircraft will be fielded to both corps 
and division cavalry and attack battalions/squadrons with the assigned missions of armed 
reconnaissance and attack, both with an embedded air to air capability. Mission profiles 
include day, night, reduced visibility, from ships, and in desert and non-desert 
environments. The aircraft is shown in Figure 1.1. 

Comanche will contribute to combat effectiveness by its ability to fight for 
information and disrupt the enemy's tactical intelligence collection effort in the counter 
reconnaissance role. Key attributes of the aircraft include a five-blade bearingless main 
rotor , FANTAIL anti-torque system, composite materials, Low Observables (LO), an 
advanced pilotage and targeting system, and an integrated digital mission equipment 
package (MEP). These digital systems will allow it to serve as a forward data fusion 
center and provide near real time information to commanders at all levels. 

Designed for two crew members, it is capable of being operated by a single crew 
member from either station. During missions, the pilot-vehicle-interface (PVI) is designed 
to support low workload, high situational awareness, human error resistance, and 


maximum mission productivity. [Ref. 1] 





Figure 1.1: US Army RAH-66 Comanche Helicopter. 


The Boeing Defense and Space Group's Helicopter Division of Philadelphia, 
Pennsylvania and United Technologies’ Sikorsky Aircraft Division of Stratford, 
Connecticut were awarded the demonstration/validation (DEM/Val) phase of the contract 
for the Comanche. These contractors divided the aircraft into two sections for design and 
fabrication. Boeing has responsibility for the tail section of the aircraft. Sikorsky has 
responsibility for the forward portion of the aircraft fuselage, to include the integration of 
both sections. 

Currently, the first prototype is undergoing extensive flight testing at the Sikorsky 
test flight facility in West Palm Beach, Florida. Focused on "envelope expansion", the 
aircraft had flown approximately 91 hours, reached 171 knots in forward, level flight, 
achieved bank angles of 45°, and sideward flight of 70 knots as of May, 1998. Prototype 


#2 1s scheduled to begin flight testing in 1999. 


5. PROBLEM DEFINITION 


Durning flight testing the empennage of the aircraft experienced buffeting where the 
horizontal and vertical tail vibrated excessively in certain flight regimes. Observable from 
both the chase aircraft and on flight data records for strain gauges and accelerometers, 
these high buffet loads created higher than anticipated fitting loads, particularly in the tail 
region, and resulted in vibrations in the crew stations and at the nose cone where the 
targeting sensors are located. Significant effort has been devoted to identifying the 
sources of excitation and the nature of the structural response. One such method is shown 
in Figure 1.2. Here, strips of maternal were placed on the empennage to allow for better 


visualization of both the airflow and vibration of the structure. 





Figure 1.2: Comanche Prototype Test Flight. 


Boeing and Sikorsky have also conducted several wind tunnel experiments to 
better quantify the flow coming from the aircraft fuselage. The first was completed in 
April 1997 at NASA Langley. Two follow-on tests were executed at the United 
Technologies’ Research Center in Hartford, Connecticut and then Boeing Defense and 
Space Group’s Helicopter Division in Philadelphia, Pennsylvania. These tests consisted of 
changing aerodynamic surfaces to modify the flow impinging on the tail from the fuselage. 
Additionally, the flight test aircraft had several structural modifications to stiffen the tail 
section; notably a bracket to secure the horizontal and vertical tail, use of a ngid flap on 
the right trailing edge of the horizontal stabilizer, and the adding of tip weights to the 
horizontal stabilizer. 

This research determines the location and magnitude of empennage loads caused 
by the random excitation. Because the nature of the excitation is assumed to be a random 
process, spectral analysis is used. To obtain the loads, a three step process 1s utilized. 
First, from aircraft differential pressure transducers and accelerometers, the spectral 
content of the excitation and response is determined. Then, using a NASTRAN model 
modified to replicate the flight test aircraft, frequency response functions are obtained 
between selected points on the aircraft's tail and the accelerometers. Finally, using this 
information, estimates of the airloads are calculated. Having provided information on the 
nature of the excitation, structural modifications can then be explored to reduce the 
response of the airframe to the excitation. Potential methods for accomplishing this may 
be based on shifting airframe natural frequencies, increasing structural damping, or adding 


vibration control devices. 


Il. HELICOPTER VIBRATIONS 


A. GENERAL 


As with all helicopters, the RAH-66 experiences vibrations. These forces can 
significantly impact crew effectiveness, imaging and targeting performance and airframe 
and equipment reliability and maintainability. With the use of complex digital avionics 
systems, there is a greater emphasis on vibration reduction. As such, the analysis of 
dynamic loads and modification of structures to create a low vibration environment is 
essential. The first step in reducing vibrations is to identify the source of excitation. Once 
that is known, the aircraft can be modified so that the structural frequencies are offset 
from the driving frequencies. 

Helicopter vibrations fall into three categories. First, there are vibrations due to 
rotor excitation, which are at frequencies that are integral multiples of the rotor’s 
rotational speed. Second, vibrations due to random aerodynamic excitation where the 
frequency being observed is a natural frequency of the structure being excited. And 
finally, self-excited vibrations, such as flutter and ground resonance. [Ref. 2] 


t Vibrations at Integral Multiples of Rotor Speed 


The primary source of helicopter fuselage vibrations is the rotor system. Here, 
harmonics of aerodynamic loads on the blade give nse to vibratory response of the blade. 
Since the blade is restrained at the root, the blade responses result in root shears, which 
feed from the rotor head into the fuselage as vibratory shears and moments. As the forces 


go from the rotating to the fixed fuselage system, the rotor system tends to act as a filter. 


For an N blades rotor system, the frequencies which will filter through are those at N and 
2N/rev. Since the amplitudes of the lower harmonics of blade loading are greater than the 
higher, the N/rev vibration of the fuselage is most critical. It can be shown that N/rev 
fuselage vibrations in the fixed system are the result of the N-1, N, and N+1/rev vibratory 
response of the blades in the rotating system. 

In addition to the N/rev or NP vibration, there are secondary excitation frequencies 
which also must be taken into account. Any unbalance of the rotor system, such as blades 
out of track, will give rise to 1/rev or 1P excitation of the fuselage. Dissymmetry of blade 
dampers may cause troublesome 2/rev excitation. Also important is the 1 x tail rotor 
rotational speed. 

Forward speed plays a critical role in vibrations. Figure 2.1 is a typical plot of the 
variation in vibration levels as a function of airspeed. The increase in vibrations levels in 
the 30 to 45 knot range can be attributed to blade vortex interaction due to the proximity 


of the rotor to its wake and the resulting non-uniform induced velocity field. 
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Figure 2.1: Typical Helicopter Vibration Levels Versus Airspeed. After Ref.[2]. 


In the high speed flight regime, the increase in vibration with increase in airspeed is 
attributed to two sources. First, the aerodynamic excitation due to impingement of the 
rotor’s wake on the helicopter’s horizontal empennage. Second, amplification of blade 
section aerodynamic loads associated with increased cyclic pitch accompanied by 
increasing differences in dynamic pressure per blade revolution. Figure 2.2 shows the 


non-distorted rotor wake geometry for both low and high forward flight speeds. 


LOW SPEED ROTOR WAKE HIGH SPEED ROTOR WAKE 





Figure 2.2: Rotor Wake Geometry. From Ref.[2]. 


2. Random Aerodynamic Excitation 

This type of vibration is often associated with the turbulent downwash from the 
rotor, which impinges on the fuselage and its components. Similar to white noise, this 
excitation contains many harmonics that can excite fundamental fuselage or empennage 


modes. Often, an interesting characteristic of this type of excitation is that it 1s present 


only within a certain range of airspeeds. Successful solutions of this problem have 
included modifying the structure excited and by altering the flow from the rotor. 


3. Self-Excited Vibrations 


A self-excited vibration is one that results in divergent oscillations where the 
system damping is negative. In the case of negative damping, the damping force, which is 
now a driving force, does positive work on the system. The work done by this force is 
converted into the additional kinetic energy of the increased vibration. A self-excited 
vibration cannot exist without an external source of energy. In the case of flutter, the 


airflow provides the source of energy. 


B. VIBRATION SPECIFICATIONS 


Vibration specifications for Army aircraft are governed by Aeronautical Design 
Standard (ADS) -— 27 and MIL-STD-810D [Ref. 1]. Human factors vibration 
requirements are expressed by ADS-27 in terms of an intrusion index which is the square 
root of a weighted sum of squares of four vibration components in each direction. The 
purpose of this is to bring together the twelve most significant crew vibration components 
into a single number for evaluation. A separate specification exists to capture the lower 
frequency one per main rotor revolution, or 1P, vibration. 

The changes in vibrations that result for airspeed and maneuvering are addressed in 
ADS-27 by defining limits in terms of flight regions. There are three defined regions, with 
the first being*normal flight up to cruise speed, then maneuvers having a duration less than 


3 seconds, and lastly maneuvers having a duration greater than 3 seconds. There is also a 


region for gunfire vibration. The specifications for the RAH-66 crew seat are listed in the 
Performance Weapon System Specification [Ref. 1] and are shown in Table 2.1. 
Additionally, Figures 2.3 and 2.4 display the currently measured intrusion index in the 


Pilot’s and Copilot’s crew station versus the specification. 


15 
40 


Flight Region III 


Weapons Firing Increment 


Table 2.1: Crew Seat Intrusion Index and MR1P Vibration Limits. 












Because of the frequency bands that are used to calculate the intrusion index, any 
random, broad band, excitation occurring between main rotor 1P and NP 1s not totally 
captured. Therefore, the intrusion index does not indicate the contribution of the broad 


band excitation to total response. 
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Figure 2.4: Intrusion Index — Copilot’s Station. 
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C. COMANCHE VIBRATION SOURCES 
12 Main Rotor Hub Forces and Moments 


As noted earlier, the primary source of helicopter fuselage vibrations is the rotor 
system. As the forces and moments are transmitted from the rotating to the fixed fuselage 
system, the rotor system tends to act asa filter. It can be shown for an “N” bladed rotor 
system that, assuming all blades see the same loading at the same azimuth position, only 
NQ frequency will be seen in the fixed system. This NQ frequency in the fixed system can 
arise from either N-1, N, or N+] excitation in the rotating system. Figure 2.5 is a 


schematic of the hub showing the mechanism of blade-fuselage coupling. 





Figure 2.5: Blade-Fuselage Coupling. From Ref.[2]. 


1} 


With the coordinate system given in Figure 2.5, forces in the rotating blade 


coordinate system can be related to forces in the fixed fuselage using the following 


transformation: 
S 
O -—j 0 0 0 oe J 
S. 2 De, 
S : 0 2 0 0 0 2 S £04 
2 2 s 
5S; |_|] 0 0 0 0 | >a 
WS ae 0 7 =e) 0 | Scos 
M. = CO Ce 0 v2 


Of {S; }= ARs, } 


where [A] is the transformation relating forces in the rotating n, € system to forces in the 
fixed X, Y, Z system. 


Thus for the five-bladed Comanche helicopter: 


e 4/rev and 6/rev flapwise blade root shears result in 5/rev hub pitching and 
rolling moments in the airframe. 


e 5/rev flapwise blade root shears feed into the fuselage directly as 5/rev vertical 
forces. 


e 4/rev and 6/rev chordwise root shears produce 5/rev airframe hub forces in the 
fore and aft and lateral directions. 


e 5/rev chordwise root shears results in 5/rev hub yawing moments. 


2 


For the RAH-66, the nominal rotor speed is 354.9 rpm. This places the 1/rev or 
1P excitation at 5.915 Hz, with the 5P blade passage frequency occurring at 29.575 Hz. 
[Ref. 3] 

2. Aerodynamic Excitation 

Aerodynamic excitation to the airframe and empennage arises from several 


sources: 


e Main rotor wake impingement. 

e Main and tail rotor downwash (induced flow) impingement. 

e Flow from around the main rotor hub and pylon assembly structure 
e Flow from attached aerodynamic surfaces. 


e Flow from protruding irregular fuselage structure, ranging from FABS 
(forward avionics bays) to as small as a canopy or door handles. 


The interaction of this flow is a complex subject, very difficult to predict, and the 
subject of numerous ongoing research projects. Rotorcraft can have problems where 
rotor airloads affect airframe loads and, in turn, the airframe effects rotor airloads and 
blade dynamic response. Some examples: The Sikorsky H-3 had a lateral tail shake that 
had its onset at about 120 knots. The problem was solved by installing a fairing or 
“beanie” on top of the main rotor hub to smooth out the flow from this region which was 
impinging on the aft fuselage. 

In development of the Boeing YUH-61 and Sikorsky YUH-60 helicopters, both 


companies obtained a significant improvement in N/rev vibration levels by “raising the 
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mast”, that is, increasing the separation between the main rotor hub and the fuselage. In 
development of the Hughes YAH-64, the initial T-tail design was scrapped and a movable 
horizontal stabilator installed to eliminate an excessive vibration problem associated with 
the T-tail. A visual summation of the typical sources of tail excitation is shown in Figure 


2.6. 


Main Rotor Vortices 


Random Excitation 
from Flow Separation 





Figure 2.6: Typical Sources of Tail Excitation 


An experimental investigation of the interactions between a rotor and a T-tail 
empennage was conducted by Lieshman and Moedersheim [Ref. 4] at the University of 
Maryland Rotorcraft Center in 1997. Using a 1/6 scale four-bladed helicopter rotor with 
generic body shape, a T-tail empennage was placed down stream of the rotor hub as 
illustrated in Figure 2.7. Both the vertical and horizontal tail surfaces used the NACA 


0012 airfoil section with a constant chord of .203 meters. The horizontal tail was located 
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in two different positions, a high and low plane configuration, with zero degrees angle of 


attack. Additional parameters of the wind tunnel test are contained in Reference 4. 


Figure 2.7: Schematic of Rotor/Body/Empennage System. From [Ref. 4]. 


Through their analysis, they drew the following conclusions: 


e Flow behind the rotor was highly asymmetric. The two low pressure regions 
trailing from the edges of the rotor disk provided evidence of the wake rolling- 
up into two larger vortex bundles. 


e At low advance ratios the flow over the horizontal tail was stalled. The flow 
gradually reattached as the advance ratio was increased. Airloading was 
asymmetric primarily due to the difference in strength of the two trailing vortex 
bundles trailing from the sides of the rotor disk. In some test conditions, a 
sharp trailing edge peak pressure was observed at the trailing edge of the 
horizontal tail. This indicated the potential of a junction vortex onginating at 
the junction of the vertical and horizontal tail. 


e Measured time-varying pressures showed the flowfield near the horizontal tail 
was highly unsteady. Below advanced ratios of .15 the unsteady pressure 
response was low. As the wake transitioned on to the tail, unsteady pressures 
increased significantly. 


e Unsteady pressure responses are characterized by two events, unsteady lift 


produced by the rotor wake and disturbances caused by the convection of 
individual vortices. 


Is 


For the RAH-66, the rotor wake transitions above the T-tail between 40 and 60 
knots, as shown in Figure 2.5 [Ref. 3]. Beyond this airspeed range the actual wake is not 
impinging on the horizontal stabilizer. From Leishman’s research; however, there can be 
convected vortices associated with the passage of blades that can induce disturbances and 


oscillating loading conditions. 
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Figure 2.8: Rotor Wake Distance from Horizontal Stabilizer. After Ref. [3]. 


In addition to the rotor vortices, the wake behind the rotor hub constitutes a low- 
pressure sink that can draw the flow off of the upper fuselage and produce separation 
[Ref. 5]. The frequency content of this excitation is often determined by modeling the 
rotor hub as a circular cylinder where vortices are shed as a Karman vortex street. 
Through the Strouhal number, a non-dimensional constant, vortex shedding frequency can 


be related to freestream velocity and the diameter of a two dimensional cylinder. [Ref. 5] 
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The relationship 1s given by: 


where f; = Vortex shedding frequency (Hz) 
S = Strouhal number 
U = Freestream velocity (ft/sec) 


D = Cylinder diameter (ft) 


Determination of the Strouhal number from Schlichting [Ref. 5:p:32] requires 


interpolating from a small plot. Using the results from a wind tunnel experiment done at 


the Naval Postgraduate School on a two-dimensional cylinder by Small, Hebbar, and 


Platzer [Ref. 6], the Strouhal number and vortex shedding frequency can more easily and 


accurately be determined. In this experiment [Ref. 6], the vortex shedding frequency and 


Strouhal number were calculated for Reynolds numbers ranging from 2.7x 10° to 1.9x10°. 


From this, the data presented in Figure 2.9 was generated. The Strouhal number 


represents the slope of the line and can be used to calculate the frequency of shedding 


from a know diameter and given velocity. 
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Figure 2.9: Shedding frequency for the Quasi-two-dimensional Cylinder. After Ref. [6]. 


For the Comanche, at 100 knots (168.78 ft/sec) and with a rotor hub pylon of 
approximately .75 meters (2.46 ft), the calculated vortex shedding frequency is 11.827 Hz. 
As the airspeed and point of separation changes, the shedding frequency will be effected, 
as shown by the equation on the previous page. The area around the main rotor pylon is 
actually a complex array of surfaces needed to reduce radar signature. Separated flow 
could occur from any of these surfaces, with varying frequency content of the vortices. In 
this case, more broadband vortex shedding will occur. 

From these factors, RAH-66 response data should show energy associated with 
main rotor 5P, 29.75 Hz and if a rotor track and balance problem, main rotor 1P at 5.95 
Hz. There could also be energy associated with separated flow from the main rotor pylon, 
centered around 11-12 Hz, and some broad band excitation from the convected rotor 


vortices impinging on the horizontal tail. 
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lil. SPECTRAL ANALYSIS 


Figure 3.1 shows the time history of lateral excitation at the pilot’s heels. When 
analyzed over the full record length of 60 seconds, the data 1s classified as random. Such 
random time functions as this data set commonly exist in nature; a typical example being 
pressure gusts aircraft encountered in flight. The method used to characterize these 
functions is known as spectral analysis, and a discussion of the theory, from Thomson 


[Ref. 7], and techniques for digital signal analysis follows. 
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Figure 3.1: Time History of Lateral Acceleration. 


A. TIME AVERAGING 


Despite the irregular nature of the time history in Figure 3.1, most “random” 
functions exhibit some degree of statistical regularity and by applying averaging processes, 
we can characterize the phenomena. In observing the data, each record is referred to as a 


sample and the collection of samples is the ensemble. If the ensemble average is taken at a 


We, 


specific time, ¢,, and then at another time, ¢; + 7, and they remain the same, the process is 
said to be stationary. Similarly, if the ensemble averages are replaced by time averages, or 
averaging each sample, and the averages from different samples remains the same, the 
random process is ergodic. Key assumptions necessary for our analysis is that the random 
response phenomena exhibited by the aircraft is stationary and ergodic. As such, statistical 
properties can be derived from single time histories of sufficient length. [Ref. 7] 

To derive these properties a review of statistics is required. From the time 


histories, we can define the expected value (mean): 


E[x«]= lim [x@nar 


Too 0 


For the case of discrete variables, the expected value is given by the following equation: 


. 1X 
Alx|=lim— 2x, 
n—->o i=l 
We can then apply the averaging process to the mean square value, also referred to as the 
square of the expected value. 


Ex? (|= x? = lim Jx"a! 


T->0 0 
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Linking this back to basic statistics, with the variance defined as: 
oOo” = x? = (a 


es 
o’ = limp! © ~¥)*dt 

the mean square value (MSV) is simply the variance when the data has an expected value 
(mean) of zero. If thought of in terms of a sin wave, or multiple harmonic function, where 
the mean is removed (detrended) before analysis, then the variance and MSV are the same. 

Knowing these statistical properties of a given signal, the correlation, or similarity, 
between different signals can be determined. This is an important function in time domain 
analysis as it can help identify signals buried in noise or propagated in time through a 
structure. Additionally, it provides a transition into Fourier analysis and the frequency 
domain, where a structure’s spectral characters can better be analyzed. The correlation 
function is computed by multiplying the ordinates of two records at each time and 
determining the average value by dividing the sum by the number of products. The 
product is a maximum when the two signals are very similar or identical. For dissimilar 
records, the products will be negative and positive, so the sum will be small. Iftwo 
signals are identical with one shifted by some time, we can compute the correlation by 


taking just a single time record and multiplying it by itself, with the specified time delay 1. 


This is called the autocorrelation function, and is defined by the following formula. 


%, 
R(z) = E[x() x(t +7) = (x) x(t + 7)) = lim= [xOx(t +7 )dt 
To 1, 


Zh 


If the time delay t = 0, the equation reduces to the mean square value or variance. 


RO)" Om 
If the two records are truly dissimilar, we call their product the cross correlation. It is 


define it by: 
R,, (7) = Elx(Q y+ 7) = (xOy +7) 


rT 
2 


= lim= [x@Ove + t)dt 


T—-0 Tf, 


B. FREQUENCY DOMAIN ANALYSIS 


Having defined averaging processes for the time domain, the random function can 
be transitioned to the frequency domain [Ref. 7]. This will assist in relating the response 
data back to the sources of excitation discussed earlier. 


The time signal x(t) can be expressed in terms of its Fourier series: 


c=) Gre a 


Because the record 1s periodic it can be represented by the Fourier series: 


x(t) =Re>C,e”™ = =e" aie a) 
n=) n=] 


a2 


where C,, a complex function and C,, is its complex conjugate. It can then be proven [Ref. 


7] that the mean square value of the function is given by 


From this, it can be seen that the MSV can be computed in either the time or frequency 
domain, depending on which technique provides better analysis of the data. 

The mean square value (MSV) represents contributions from the entire frequency 
band of the signal. If only a certain frequency range is of interest, the MSV for that range 
of frequency, Af, is the power spectrum (PS). The power spectrum, P(7,), is then just the 
product of the Fourier coefficients from the frequency range of interest, f,, with the mean 


square value being the sum of the power spectrum across the entire bandwidth. [Ref. 7] 
pages: ae 
PSI = SOC, ae a) 
n=] 


If the power spectrum is divided by the frequency interval, Af, the result, defined 
as the discrete power spectral density(PSD) S(f,), gives the contribution to the power 
spectrum from an even more narrow frequency band. This can be extended for a 


continuous spectrum as the limiting case of S(/,) as Af > 0. 


PCF) _ CC, 


he a 


The mean square value can then be computed from the PSD by 
x = DSF, )AS 
n=} 
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The root mean square (RMS), synonymous with the standard deviation, can be found by 


taking the square root of the mean square value (MSV). 


RMS = Vx? 


As we are dealing with non-periodic functions over discrete time intervals, use of 
the Fourier integral, a limiting case of the Fourier series, is required. The Fourier Integral 


and Fourier transform of x(t) are defined respectively as 


xh) = [x(ne™*af X(f)= [xe at 
Using Parsevale theorem, the MSV of the discrete ene the time domain can 
be equated to the Fourier transform of the signal multiplied by its complex conjugate. If 
x(t) and x2(t) are real valued time functions and their Fourier transforms are given by 


Su and X2(f) then 


co © 


[xu Ox. Oat = [XNX2Naf 


-o 


If x;(t) = x2(t) = x() the MSV and PSD can be expressed in terms of the Fourier transform 


i 
Z 
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As with the Fourier coefficients, the MSV and RMS of the signal can be 
determined from the Fourier transforms. Additionally, the autocorrelation function, 
symmetric about the origin, can be related to the spectral density function via the Weiner- 


Khintchine equation. [Ref. 7] 


SC) ig 2f R(r) cos2aftdt 


The cross-spectral density (CSD), defined as the Fourier transform of the cross- 
correlation function: 
an | HOM cm adc 
Can also be determined by multiplying the Fourier transforms of the two signals of 


interest. 
S.(f)= lim>* “AY) 


Given all of these tools, the one of major interest in analyzing the flight test data is 
the power spectral density (PSD). This function of frequency gives an indication of the 
contribution to the total response, in our case accelerations, due to a discrete frequency 
bandwidth. Additionally, the cross-spectral density recaptures the phasing information 


lost when the Fourier Integral is multiplied by its complex conjugate. 


Zo 


C. DIGITAL SIGNAL PROCESSING 


Given this theory, digital signal processing methods have been devised that 
compute the PSD from discrete time histories. Vectors represent the time histories. In 
our case, the sampling rate is 300 samples per second and record length is 60 seconds. 
Therefore, the vector of the discrete signal has 18,000 data points. 

To convert this time signal into either a PSD or CSD, the first step is to decide 
how many points, N, are needed in the Fast Fourier Transform (FFT). For computational 
speed and accuracy, FFTs are done in power of two. The number of points also 
determines what the effective bandwidth of the processed information will be. Bandwidth, 
Af, is defined as 1/Time, and for a discrete signal, Time = N x At (time increment). The 
time increment, At, is equal to 1/sampling rate, or, in our case, 1/300 sec. Ifa 1024 
point FFT is utilized, Af= 300/1024 = .293 Hz. As the number of points in the FFT 
increases, the bandwidth decreases and you get finer frequency resolution. 

There is a penalty when working with random data to increasing the number of 
points in the FFT. Based on N, the signal is divided into ng records or windows and has 
any linear trend removed. As N increases, the number of windows, ng, decreases and there 


is less averaging occurring. This is demonstrated when computing the expected value. 
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To suppress leakage of the signal from one frequency range to the next, time windows are 
used that taper the time history data. This eliminates discontinuities at the beginning and 
end of the record to be analyzed. A commonly used window, that is used exclusively 
during this analysis, is the Hanning window. It is given by 
27 

w(t) = .5*(0- COs) 
and multiplied by the N data points of each ng segment. Use of the window does reduce 
the amplitude of the processed signal. In the case of the Hanning window, the output is 
375 times smaller. A correction is made in the normalizing of the PSD to account for 
this. 

Variance is introduced by tapering the signal. To reduce this, the ng segments are 
overlapped. In terms of the Hanning window, the overlapping brings the "peaks" of each 
cosine window closer to each other. This process, basically, takes the points that were at 
the end of one ng segment, where the window is near zero, and makes it the "main area" 
for the next. The FFTs are then computed for each tapered, overlapped segment. 


From the FFTs, the PSD can be determined using the following formulas: 


So (FFT), » 2° EUUFFT) )I* a0 


NT a ae = (norm(window))’ 


The expected value is the sum of the FFTs of each segment multiplied by its 


complex conjugate, or the magnitude of the FFT squared. Either way, it produces a real 


Zh 


valued function that is then divided by the number of segments, ng. As discussed earlier, 
and seen here, increasing N reduces the number of FFT averages in the expected value. 

The PSD is then determined by "normalizing" the expected value. This first 
involves multiplying by 2 as the PSD represents the one-sided autospectral density 
function. Then, the signal is multiplied by the time increment, At (or dividing by the 
sampling rate, from the definition of At) and divided by the norm of the window, squared. 
The reasons for this come from the definition of the FFT and by canceling several terms. 
Bendat and Piersol [Ref. 8:p.371] provide the exact proof; however, their notation 
requires dividing the signal by N and then adjusting for the attenuation caused by a 
tapering window. The norm of the window (its length in N dimensions) accounts for both 
elements. For no windowing, computationally done using a boxcar window, the norm’ is 
the number of points in the FFT, 1024. For a Hanning window, the norm’ = 1024 x .375 
= 384. As discussed earlier, this corresponds to the attenuation causing by using a 
Hanning window. Figure 3.2 and 3.3 provide a visual illustration of the procedure used to 
compute the autospectral density function (PSD). 

For the cross spectral density, the same steps are followed for each time history. 
Then the FFTs from each signal, x(t) and y(t) are multiplied by each other and the 


expected value is a complex function. E[FFT"] is normalized using the same method. 
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The time history is divided into segments. For computational speed and accuracy, each segment 
consists of 2" data points. The data is detrended to eliminate the stationary component. 


A tapering window, such as a Hanning, can be used to reduce the spectral leakage. To 
counter the variance caused by tapering, overlapping time histories are used. 


Hanning Window Hanning Window Hanaing Window 





Figure 3.2: PSD Computational Steps — Time Domain. 
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The FFT is computed of each tapered time history. The magnitude of the FFT is then squared. 
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Figure 3.3: PSD Computation — Frequency Domain. 
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The PSD or CSD output vector is of length N/2 with a frequency band from 
0:Af:(sampling rate/2). As an example, with N of 1024 points and a sampling rate of 300 
Hz, the maximum frequency, commonly referred to as the Nyquist frequency, 1s 150 Hz. 
The frequency vector begins at 0 and goes to 150 Hz in increments of 300/1024 or .293. 
For each point of the frequency vector, the PSD vector has a corresponding value. This 
value represents the energy associated with that frequency band. So, in this example, the 
first value of the PSD vector represents the energy present within the 0 to .293 Hz range. 

Units of either the PSD or CSD are (input units)*/Hz. In the case of accelerometer 
output, the PSD units are g’/Hz. With units like this, it is often difficult to grasp an 
understanding of what the information really means. The RMS 1s much easier to interpret, 
which has the same units of the input, in this case, g’s. To compute the RMS, the PSD 
vector is multiplied by the frequency increment, Af. This vector is then summed within the 
frequency range of interest. The resulting scalar is the MSV and the RMS 1s the square 
root. If, for example, the RMS between 5-25 Hz is needed, the PSD vector is multiplied 
by Af = 300/1024. Then, only the indices of the vector that represent energy between 5- 
25 Hz are summed, points 17-86 of the 512 total PSD values. The square root is the 
RMS and represents the total energy or response between 5 and 25 Hz. Comparing back 


to the time domain, this RMS value equates to the standard deviation. 


3) 





IV. FLIGHT TEST DATA ANALYSIS 


Time histories of accelerometers and pressure transducers from flight tests were 
analyzed using both 1024 and 4096 point FFT. For computational speed, 1024 point 
FFTs were utilized for final analysis. Several different tapering windows were also applied 
with no noticeable difference. The excitation frequencies are clearly evident and are an 
order of magnitude higher than the energy of the surrounding frequency bands. Therefore, 
even when no side lobe reduction techniques were utilized, the frequency content of the 
response was pronounced. To remain consistent with the methods used by Sikorsky, a 
Hanning window with 67% overlap was used. This did reduce the spectral leakage and 


variance associated with windowing and produced a more accurate PSD value. 


A. ACCELEROMETERS 


Data records at 60, 80, 100, 120, and 140 knots were obtained for the 
accelerometers listed in Table 4.1. The sampling rate for all signals was 300 samples per 
second with a 60 second record length. A complete listing of all test flight parameter is 
contained in Appendix A. 

The PSD of each accelerometer through the range of airspeeds was computed. In 
looking at the energy across the entire frequency spectrum, 0-150 Hz, there was 
significant energy from 0-30 Hz, and none of magnitude beyond that range. Figure 4.1 
through 4.3 display the PSD for nose cone lateral, pilot heel lateral, and copilot heel 


vertical. The programs used to compute the PSDs and CSDs 1s in Appendix B. 
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_ Accelerometers 


Table 4.1: Accelerometers. 


















As previously noted, primary excitation for a 5-bladed rotor occurs at 5P, which is 
29.6 Hz for the Comanche. As expected for this case, 5P response centered at 29.6 Hz 
dominates, especially at higher airspeeds. In contrast, the amplitude of the other frequency 
bands is not discernable. To provide more insight at response due to other than SP, 
Figures 4.4 through 4.18 displays the accelerometer response from 5-25 Hz. Within this 
frequency band, much greater detail in the response data can be discerned. Given on the 
next 17 pages are plots that indicate the frequencies of excitation along with the expected 


% 


increase with airspeed. 
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PSD Nose Cone - Lateral Response 
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Figure 4.1: Nose Cone — Lateral 5-30 Hertz. 
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Figure 4.2: Pilot Heel — Lateral 5-30 Hertz. 
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Figure 4.3: Copilot Heel — Vertical 5-30 Hertz. 


PSD Nose Cone - Vertical Response 
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Figure 4.4: Nose Cone - Vertical 5-25 Hertz. 
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PSD Nose Cone - Lateral Response 
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Figure 4.5: Nose Cone — Lateral 5-25 Hertz. 
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Figure 4.6: Pilot Heel — Lateral 5-25 Hertz. 
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Figure 4.7: Pilot Heel — Vertical 5-25 Hertz. 
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-4 PSD Pilot Floor - Vertical Response 
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Figure 4.9: Pilot Floor — Vertical 5-25 Hertz. 
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Figure 4.10: Copilot Heel — Vertical 5-25 Hertz. 
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Figure 4.13: Copilot Overhead — Lateral 5-25 Hertz. 
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PSD Vertical Pylon Top - Lateral Response 
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Figure 4.14: Vertical Pylon Top — Lateral 5-25 Hertz. 
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PSD Horizontal Stab Left Tip - Vertical Response 
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Figure 4.15: Horizontal Stab Left Tip — Vertical 5-25 Hertz. 
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PSD Horizontal Stab Right Tip - Vertical Response 
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Figure 4.16: Horizontal Stab Right Tip — Vertical 5-25 Hertz. 
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PSD Horizontal Stab Left Tip - Longitudinal Response 
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Figure 4.17: Horizontal Stab Left Tip — Longitudinal 5-25 Hertz. 
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Figure 4.18: Horizontal Stab Right Tip — Longitudinal 5-25 Hertz. 
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From a review of the plots, there are three frequency bands that continuously 
display energy. Energy in this context is based upon measured vibration response data as 
either acceleration in g’s or PSD in (g’s)*/Hz. The first band is centered around 5-8 Hz, 
the second from 10-13 Hz, and the last 22-25 Hz. From the PSD data, the RMS for each 
airspeed was calculated. Figures 4.19 through 4.22 show the computed RMS from the 
nose cone, pilot’s station, and copilot’s station. The RMS including and not including 5P 
is plotted along with both lateral and vertical response. The 5P vibration adds a significant 
amount to the total response and produces a sharp rise in accelerometer readings 
acceleration at airspeeds above 100 knots. The total response for the frequency band of 
5-25 Hz begins to level off at 100 knots, the slope being much less than the RMS with 5P 
included. This excitation, not including 5P, is less prevalent at lower airspeeds and 
significantly increases from 80-100 knots. Past 100 knots, this frequency band energy 
begins to level off while 5P vibration levels increase significantly. By plotting both lateral 
and vertical RMS, comparisons as to the more dominant response can be made. Lateral to 
vertical ratios typically range of the order from 1.25:1 to 2.0:1. In all cases except for the 
pilot’s floor, lateral response dominates. To capture the maximum excitation outside of 
the 5P range, and minimize 5P response, a more detailed analysis is presented of the PSD 
at 100 knots. 

Figures 4.23 through 4.28 present the PSD for 100 knots. Again, both lateral and 
vertical responses are plotted together. As was seen in the RMS calculations, the lateral 
response is dominant. Additionally, the three frequency regions identified earlier are 


evident. 
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Figure 4.19: Root Mean Square — Nose Cone. 
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Figure 4.20: Root Mean Square — Pilot Heels. 
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Figure 4.21: Root Mean Square: Pilot Floor. 
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Figure 4.22: Root Mean Square — Copilot Station. 
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PSD Nose Cone - 1024 Point FFT 
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Figure 4.23: Nose Cone — 100 Knots. 
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Figure 4.24: Pilot Heels — 100 Knots. 
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PSD Pilot Floor - 1024 Point FFT 
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Figure 4.25: Pilot Floor — 100 Knots. 
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Figure 4.26: Copilot Station — 100 Knots. 
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Figure 4.27: Horizontal Stab Left Tip — 100 Knots. 
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Figure 4.28: Horizontal Stab Right Tip — 100 Knots. 


58 


To quantify the amount of response due to these three frequency bands, or bins, 
RMS values for the discrete ranges of 5-8 Hz, 10-13 Hz, and 22-25 Hz were calculated. 
Each RMS value was summed and compared to the total RMS, not including 5P response. 
The data is presented in Table 4.2. Although these three frequency bins represent only 9 
Hz, or 45% of the total bandwidth, they constitute from 50% of the total response to as 
much as 90%. Comparing the bins, the more dominant frequency range was 10-13 Hz, 


with 60% of the accelerometers having their peak response in this region. 













_____| _Frequeney Bins __ % Total 


10-13 Hz| 22-25 Hz|_RMS 
0.015 


ota 

0.8355_ 

07390 J osser_| 61.18 — 
Table 4.2: RMS of Frequency Bins — 100 Knots. 
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5: FLIGHT TEST PRESSURE MEASUREMENTS 


The differential pressure transducer data records that correspond to the 
accelerometer data are listed in Table 4.3. The sampling rate remained 300 samples per 


second with a 60 second record length. 


- Differential Pressures 


__ | Differential Press Pressure, ure. 7% c chord, 30% 
up Vertical Tail 


Differential Pressure, 7% chord, 60% 


VERTDP2 


ee up Vertical Tail 


Differential Pressure, 7% chord, 60% 
Span Left Honzontal Tail 
Differential Pressure, 7% chord, 60% 
ee Span Right Horizontal Tail 

Table 4.3: Differential Pressure Transducers. 


LHORZDP3 





The frequency content of the differential pressures was determined and compared 
to the response PSDs, findings of Leishman and Moedersheim [Ref. 4], and the 
approximated shedding frequency of the rotor hub pylon. The analysis also provided 
insight into which areas of the structure are excited. This aided in focussing the 
NASTRAN analysis, deriving the transfer functions between the random aerodynamic 
loading and the response locations, and correlation with mode shapes. This will be 
addressed in Chapter VI. Figure 4.29 shows the variation of differential pressure from the 


left horizontal stabilizer as a function of airspeed from 60 to 140 knots. 
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Figure 4.29: Differential Pressure - Left Horizontal Stabilizer. 

As with the response data, the magnitude of the 5P differential pressure (29.85 Hz) 
is an order of magnitudes greater than excitation in any other frequency band. To gaina 
better appreciation for the vibratory pressures at other than 5P, the transducer output is 
plotted for the 5-25 Hz range in Figures 4.30 through 4.33. 

Given are three dimensional plots of airspeed, frequency and magnitude of the 
PSD. The asymmetry between the left and right horizontal tail is expected due to the 
angle of the vertical tail. The right horizontal tail exhibited pressure spikes at the 


previously noted discrete frequency bands of 5-8 Hz, 10-13 Hz, and 22-25 Hz. 
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Figure 4.30: Differential Pressure - Right Horizontal Tail. 
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Figure 4.31: Differential Pressure - Left Horizontal Tail. 
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Figure 4.32: Differential Pressure — Vertical Tail. 
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Figure 4.33: Differential Pressure - Vertical Tail. 
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Unlike the discrete excitation of the nght horizontal tail, the left honzontal tail has 
broad band excitation from 11-19 Hz. Both span locations on the vertical tail show 
symmetry in the shape of broad band excitation, similar to that of the left horizontal. 
There is a decrease in magnitude, on the order of 12:1, as you move up the vertical tail, 
from the 30% to 60% span location. 

Figures 4.34 through 4.35 are the PSDs of differential pressures at 100 knots. 
Both horizontal and vertical transducers are contained in the same plot for better 
comparison. For the horizontal tail, the magnitude of the pressures seen on the left is 
much great than that of the right. The energy is peaked within the 11.5 to 14 Hz band, 
with an additional peak at 18 Hz. The right displays its peak pressure at 6 Hz, which 
corresponds with the 1P of the aircraft. 

The vertical tail, at both the 30% and 60% span location, have the same 
frequency distribution as the left horizontal, 11.5 to 14 Hz. As seen earlier, the magnitude 
of the pressure continues to diminish as you move up the vertical tail. Although the 
pressure distribution between the vertical tail and left horizontal is similar, the pressures 
seen at 30% span of the vertical are 3.5 times greater than that of the left horizontal. The 


pressures at the 60% span of the vertical and the left horizontal are of equal magnitude. 
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Figure 4.34: Differential Pressure Horizontal Tail - 100 Knots. 
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Figure 4.35: Differential Pressure Vertical Tail - 100 Knots. 
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Linking the pressure distributions back to response data, RMS values were 
calculated for the discrete frequency bins used earlier. The unit for the RMS of the 
pressures is psi. Table 4.4 presents both the magnitudes of pressures as you move up the 


vertical tail and along the horizontal as well as the dominant frequency band. 






Root Mean Square - 100 Knots - psi 


re a 
7 Ch 30% Sp Vert ~ 0.0503 


Table 4.4: RMS of Differential Pressures — 100 Knots. 












Understanding the flow over the empennage, even for the limited amount of 
transducers, assisted in focussing the analysis of the structural response. The following 
chapter presents the structural analysis and correlation between aerodynamic forces acting 


on the surface and response at accelerometers. 
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V. STRUCTURAL DYNAMICS 


The dynamics of the multiple degree of freedom (MDOF) Comanche structure 
govern the transfer of random loads to fuselage response. In order to understand the 
MDOF system, the single degree of freedom case, also called a single modal oscillator, 
must first be studied. From these solution methods, the dynamics of the complex structure 
can then be quantified using a systematic modal analysis approach. Of particular interest is 
the response of the system to harmonic, or frequency dependent excitation. The quantity 
that describes this is called the frequency response function (FRF) or, from the Laplace 
domain, the transfer function evaluated at (j@). Working from a vibrations perspective, 
the FRF for a single DOF system will be computed and then applied to a MDOF 


Comanche structure. 


A. SINGLE DEGREE OF FREEDOM (SDOF) OSCILLATORS 


The response of a single degree of freedom system to a forced harmonic excitation 
can be determined from the equation of motion that describes the system. The SDOF 


oscillator is characterized by 


mx +. cx + kx = F cos(Q?) 
where 
Total (modal) mass participating in mode of vibration 
Total (modal)stiffness participating in mode of vibration 
Total (modal) damping participating in mode of vibration 
Magnitude of forcing function 


Oo eS 
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The equation can then be divided by m and is in the generalized form, 


- ae F 
¥+2x+"x= —cos(Q7) or ¥+260,X+@;x = SO 
mm m m 


The total solution is comprised of two parts: 
x(t) = X homo (Z) am X part (f) 


The homogeneous solution, also known as the transient or free response is 


Xnomo (C1) = e °°" (C, cos(w,t) + C, sin(@,t)) 
where the constants C; and C2 are determined by applying initial conditions. The 
particular solution is found for a specified forcing function. In the case of a cosine forcing 


function, the solution is 


Xan (2) = Acos(Q4¢t —- PD) 


and the system responds at the frequency of excitation, 22. 
To solve for the amplitude and phase of the response, A and @, complex algebra is 


used and the general forcing function becomes 


| x 


F(t) = Fe = F(cos(Q¢) + jsin(Qz)) where F = 
m 


With this definition, the general equation for the forced response becomes 
¥+26@,%+ 02x = Fe™ 


a 


and if the assumed solution for the particular solution is 


vaa(e) = Ae 
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the necessary derivatives can be taken to find A in terms of known quantities 


(-Q? +26QQo, j+o2)A=F 


dividing by the natural frequency squared, @,” and grouping the real and imaginary parts, 


(1-226 2s} Z(Q)e¢A=F 
@ @ 


n n 


with the force amplitude redefined as 


ofa; 
@ k 


ni 


The terms inside the brackets represent the complex frequency dependent constant, 
Z(Q), known as the impedence, which is a function of the frequency Q. The equation can 


also be wnitten in the form 


Converting to polar form and dividing magnitudes and subtracting angles, a solution for A 


IZ seas (Q) + JZ snag (QA = F 


is obtained 


ja gi) 
Z seu (2)? + Zipag (OY? 
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The modal parameters defined in our equation of motion can be substituted in to yield 


Ae a ee e? 


rr 
QO, Q,, 


The particular solution is then defined as 


where @® = tan” 


Xoon (t) = [Ale = |Al((cos(Qr - ®) + j sin(Qr - ©)) 
As the forcing function was a cosine function, only the real part is taken and the 


general solution 1s 


X sant (ft) = |A\(cos(Qr — ®) 
The amplitude of the steady state response is given by 


[4| = 


If the amplitude is then normalized by the static displacement F/k, the frequency response 


function (FRF) 1s found 


1 


Cae 
,, QO, 


H(Q) = 


The FRF can then be used to relate any forcing function at a given frequency to the 
output. If the forcing function is defined as F(©) and the output, say displacement, is 


Y(Q), then the relationship between the two is given by 


¥(Q) = H(Q)F(Q) 


B. MODAL ANALYSIS 


Knowing how to solve for the FRF of a SDOF system, these methods can be 
applied to the MDOF systems used to represent complex structures. In order to solve the 
coupled set of equations associated with MDOF systems, they first must be decoupled into 
single DOF modal oscillators. Each modal oscillator is then solved in the procedure just 
described and then combined to form the total solution. This method 1s called mode 
superposition and is consistent with the matrix methods currently used for most structural 
analysis. [Ref. 9] 

If the MDOF structure is defined as a set of NxN coupled equations represented 


by 


Mt} +[C]}+ [Khix} = {F()} 


The solution to this equation is first found by solving the homogeneous eigenvalue 
problem with no damping, C. Substituting an assumed solution of the form {x}={o}Qe’” 


and dividing by the mass, the eigenvalue problem 1s defined as 


TI 


My" [k]- o? [1] ko} = (0} 


The non-trivial solution of the eigenvalue problem produces a vector of length N 
that contains the natural frequencies squared, @,” with units of (rad/sec)”. The mode 


shapes, or eigenvectors {}, can be easily solved for and combined into the modal matrix. 
[o]= [fo fiw? fo? j-10* J 


The modal matrix, although not an orthogonal matrix, is orthogonal with respect 
to the mass and stiffness matrix. As such, the mass and stiffness matrices can be 
transformed by the modal matrix to diagonal matrices that decouple the equations of 
motion and yield N - SDOF modal oscillators. The coordinate system must also be 


transformed to modal coordinates, g 


ix} = [Oka 


Pre-multiplying by the transpose of the modal matrix and substituting for x, the equation 


of motion, minus the damping becomes 


[>} IMJoha}+ [oP K]okaj=lol FS 
Defining the modal mass, stiffness, and force as 


* fol MJo]=|s]; [of [kJo]= [RK]; F f= fof 


oe 


The modal equation of motion is 
IMfa}+ [Kos = 3 


These diagonal matrices now represent the N modal oscillator SDOF equations of motion. 


A solution for each can easily be obtained. 


Damping is re-introduced on a modal basis. As each differential equation is 
solved, the damping associated with that particular natural frequency, or mode, is applied. 
This is consistent with the techniques used in parameter identification where damping 
ratios are determined on a modal basis. The response in physical coordinates is found by 


multiply each modal solution by the modal matrix 


{x}=[®]tg} 


If interested in the frequency response function (FRF) between only certain points, 


i andj, on an MDOF system, the following method is used 


p" modes oo’ 
Hi; (02) = »; > 
r=] Qo, Boe, + J 2Gaa 
Once the FRF is known for those points, the magnitude of the response, A, due to 
a harmonic excitation of amplitude F and at a driving frequency of 2, can be determined 


by multiplying the FRF and F. 
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C. RIGID BODY MODES 


Unconstrained systems have rigid body modes. This occurs when the structure 
translates or rotates without any deformation. Numerically, this appears in the solution as 
natural frequencies, ,, of zero. The reason for this is based on the relationship between 
the stiffness matrix, K, and the strain energy of the system. If the MDOF system moves 
without deforming, no strain energy 1s produced. The stiffness matrix, K, becomes 
positive semi-definite with a determinate of zero. As such, the solution to the eignevalue 
problem has as many trivial solutions, zero, as there are rigid body modes. For an aircraft, 
there are six rigid body modes, which is the case for the dynamics of any body in three- 


dimensional space. 


D. FINITE ELEMENT MODELING 


To determine the frequency response functions using modal analysis, the 
Comanche structure was modeled using the software package NASTRAN. The National 
Aeronautics and Space Administration (NASA) initially development NASTRAN in the 
1960s. NASTRAN is an acronym for NASA STRuctural ANalysis. Originally written in 
FORTRAN, it uses the finite element method (FEM) to analyze structural models [Ref. 
10]. The foundation of this is the discretization of the structures stiffness and mass. The 
FEM provides the basis for algorithms that can efficiently analyze complex structures such 
as the Comanche. A detailed discussion of the Finite Element method and its relationship 
to structural ‘dynamics is contained in Reference 11. NASTRAN version 70 was used for 


this analysis. 
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The results generated using NASTRAN were viewed using PATRAN 7.0, an 
integrated computer-aided engineering (CAE) tool. The menu-driven graphical interface 
makes model manipulation easy and aids in visualization of mode shapes and model 


geometry. 
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VI. RAH-66 NASTRAN ANALYSIS 


A NASTRAN finite element model of the RAH-66 was provided by Sikorsky. This 
model, yrahj157_reinf_2kg.dat, was modified by engineers at Boeing to replicate the 
configuration of the prototype aircraft. The name of the model is very important as 
Sikorsky manages many configurations and models, each with different stiffness and mass 
properties. The model used in this analysis included additional stiffening to the tail of the 
aircraft and 2.2 kg tip weights on both sides of the horizontal tail. With this, the natural 
frequencies, mode shapes, and frequency response functions (FRF) could be determined 
that most closely represented the actual aircraft. This was critical as these FRFs were later 
used with accelerometer and pressure transducer data to compute estimates of the loads. 


A PATRAN representation of the NASTRAN model is shown in Figure 6.1. 


A. NATURAL FREQUENCIES 


Hover frequency response testing of the test flight aircraft was conducted at West 
Palm Beach to support aeroservoelastic analysis. The purpose was to identify actual 
natural frequencies and modal damping. Through flight control frequency sweeps from 5- 
15 Hz, actual fuselage natural frequencies and percent critical damping were determined 


for selected modes. 


a 





Figure 6.1: RAH-66 NASTRAN Model. 


Rather than adjusting the stiffness and mass properties of the NASTRAN model to 
return the actual natural frequencies, engineers at Boeing demonstrated that if the mode 
shapes from the original NASTRAN model were used along with adjusted natural 
frequencies, the resulting displacements and accelerations matched actual test flight 
response data. The method used to accomplish that was to first ran the NASTRAN 
analysis deck. Then, during a restart using DMAP commands, the eigenvalue matrix 
(representing NASTRAN natural frequencies) was altered to actual frequencies. The 
resulting accelerations from application of harmonic loads closely matched actual response 
due to harmonic sweeps of the tail fan. Table 6.1 displays the NASTRAN natural 
frequencies and actual natural frequencies and damping, for selected modes, as determined 


from flight test. 
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Freq (Hz) | 2% Critical 

eee secure ee 
| 8 | (7.77 ___| Fuselage Vertical Bending | __7.70_ | 2.34% 
| 9 | 853  |PilotVert/Pitch =| | 
Pale Sewer nlomealvaw | ae 
Pa oss espilot Vewritch | eee 
Ce CS copilot lay @ ao al | 
| 15 | 13.67 __| Tail Rotor Drive ShaftMode | | 
2™ Fuselage Vertical ae 
Dit 4 ieuilsut 
| 18 | (15.75 | TailRotorDrive Shaft | | 
| 19 | 15.95 | FuselageHead Vertical =| | 
| 20 | 1607 |2™FuselageLateral | | 
| 21 | 16.55 | PilotSeatLateralRoll | | CC 
| 22 | ~=16.84 ~~ | MainXmsnLat/Roll | 
| 23 | = 17.10 | TailRotorDrive Shaft || 
| 24 | ~—=:17.82 | TailRotorDriveShaft | | 
| 25 | 18.24 =| AmmoBay TC 
| 26 | 18.82 ~~ | CopilotLat/Rol | 
| 27 | 1889 ~~ |MainXmsnLat/Roll | 
| 28 | 2045  ~|FuselageMode | 
eee ee Ponce |) eee 
| 30 | 21.56 | StabVerticalBending || 
31 | 21.78 | TailRotorDriveshat || 
| 32 | 2180 |HeadLateral | CC 
| 33 | 21.87 | TailRotorDriveshat || 
34 | 23.05 | CopilotSidearm || 
| 35 | 2402 | WeaponsDoor | iL Cis 
| 36 | 2404 ~|Fuslagge . = | CT CidC 
| 37 | 24.23 | WeaponsDoor || 
| 38 | 2453, |[BCUMode = Cd 
| 39 | 2461 | PilotSeat/Sidearm | CTC 
| 40 | 2474 [Canopy CT CT 


Table 6.1: Aircraft Natural Frequencies. 


Ty 


B: MODE SHAPES 


Selected mode shapes for free response from the original NASTRAN solution are 
contained in Appendix C. Again, as the flight test derived natural frequencies were only 
slightly different from the NASTRAN frequencies, it is assumed that the NASTRAN 
mode shapes are accurate. Three critical mode shapes are shown in Figure 6.2 through 
6.5. The first two mode shapes were identified in test flight and have large lateral 
response in frequency bins previously discussed. The third mode shape is outside of the 
frequency range; however it correlates well with the response seen in the 22-25 Hz range. 

Figure 6.2 is the aft fuselage lateral torsion mode. This motion is consistent with 
the accelerometer output in the 5-8 Hz bandwidth. Referring to Table 4.2, correlation can 
be made between the flight test output and the free response. There is large response in 
this frequency bin from the accelerometer at the top of the vertical stabilizer. This same 
motion occurs in the 6 Hz mode shape. 

The horizontal stabilizer yaw mode at 11.6 Hz, displayed in Figure 6.3, has 
significant lateral response at the base of the vertical tail and at both the pilot station and 
nose cone. The RMS of accelerometer output from both the pilot station and nose cone 
have the largest magnitude in the 10-13 Hz bandwidth. 

There is no mode shape within the 22-25 Hz bandwidth that correlates well with 
the flight test data. Figure 6.4 shows the horizontal stabilizer bending mode at 26.02 Hz. 
In the baseline NASTRAN model, yrahj170, the frequency is 22.08 Hz. This response is 
consistent with flight data and is characterized by vertical motion of the horizontal 


stabilizer and lateral response of the CPG station and pilot overhead. 


80 


s 
Ne 


Silsoo). ts a te ya: SNA WAL cle ARLULLWUVELERVOUEDEDEDRRR RRR RDR LRRD RR UR RDB BORER RRR 


aTly 


Set" 
O09 | 
PBL 
L706 








SEE POO 





"c662 
968 
tube 


| L67 
A sort 


1686S 











ZH8H O09 ‘Spo, [euoIsIO] jesoye7] ebejasny yy 





Aft Fuselage Lateral Torsion Mode Shape. 
81 


. 
. 


Figure 6.2 


oT] e988 
AQ 992 


Orb 
e25 | 
80Gb 
O62 z 
9% 
280k 
vere 
GLa 
OBI > 
LEG 


bie 
* 


RG ores 


ESaeOD 


e's 


BER he RDN 


Loe 
goar 
Bh 
82} 


1606" 
4 U6e2 


192 


"| zoe 


PES E 


ey SIRE 
~ 1 O6Ey 


PO PPDIPAS 


LES 
BGG’ 
OPES 
talG% 








MED be SMT OR CT SSS 





ZU8H O91} ‘SPO MBA QeIS |B]UOZOH 


| Stabilizer Yaw Mode Shape. 


nzonta 


Ho 


Figure 6.3 


82 





Wes ASS E KY 


“SBT 
“e] 61s 
SS 201 
PES "L 
1d g 
QxG'2 
4 9606 
~~ Be 
“1 O40 
ay BiS's 
P{seos 
P| B63'S 
o01's 


ZA £099 
RG pire 


wed L 


ARRRRREVERVEDVURVRVVRVARV RUG RR DR VURU RUUD RR ULLAL RRR U RRS RRR US RUN URR UU R RRR R RRR R RUN R RRR RNR RS WRUAVDDDVERARERREGLRRD REAR RAE DDT DENT 


> 
bl) 
A 


ae 


Pe 


‘3 
oo 
& 
>< 


s 
: 

. 

a ae 
. 

= 

wee 


, . . J 
Ot SAN@rt 





Sa “ EES A 
peemconcomaaa ys) Derr se oie ial arate Sse Wipe i= Bh 9-1 
wr 2 
Sh ne r aoe, Fe 


ZU8H Z0'92 ‘epowy Bulpuag [eolue, ges jeIUOZLOH 








Vertical Bending Mode. 


1Zer 


1] Stabil 


Horizonta 


Figure 6.4 


83 


c FREQUENCY RESPONSE FUNCTIONS 


Frequency response functions were determined between the horizontal and vertical 
tail and accelerometers. To compute the FRFs, accelerometer locations had to be 
correlated to NASTRAN enids. Sikorsky provided a list, contained in Appendix D, of 
accelerometer locations and suggested NASTRAN grids. A partial list of the NASTRAN 


grid points used 1s shown in Table 6.2. 


Table 6.2: Accelerometer Locations in NASTRAN Model. 





A one pound force was applied at these nodes with a frequency sweep of 5-25 Hz. 
Accelerations, in g’s and in all three directions, were found at each accelerometer location. 
The FRFs were wnitten to a .pch file, converted into matrices compatible with MATLAB, 


and stored as .mat files for use in solving for PSDs and CSDs of loads. 
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a" Horizontal Tail 


The NASTRAN nodes along the quarter chord of the horizontal tail are shown in 
Figure 6.5. Only one loading condition, loads applied in the z direction, was used. In 
addition to determining the FRF at each grid point along the quarter chord and the 
accelerometer locations, a chordwise investigation was conducted. At five different span 
locations, the FRF was found for each node along the chord. This provided no additional 
information as the FRF varied minimally from leading edge to trailing edge. Therefore, it 


is not included. 


Horizontal Tail Loading Condition 


Vertical Force Left Tip C/4 





Figure 6.5: Horizontal Tail Nodes and Loading Condition. 
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Figures 6.6 and 6.7 show the FRF between the horizontal stabilizer and the pilot 
heels lateral and vertical. Although the FRF was determined from 5-25 Hz, the magnitude 
beyond 20 Hz was extremely small. For better resolution in the plot, the axis was reduced 
to 5-20 Hz. The span axis represents every node along the quarter chord, with the nght 
tip defined as 0% span and the left tip representing 100%. As expected, the peak FRFs 
are aligned with the natural frequencies. 

In Figure 6.6, the lateral FRF has the greatest magnitude at 6 Hz on the right tip of 
the stabilizer. All other FRFs are less than this value. As the load was moved along the 
tail, the magnitude continued to decrease until a minimum at mid-span. This is expected 
as loads directly over the vertical tail induce very small moments, and the lateral motion 
throughout the fuselage reduces. Continuing out the left horizontal, the magnitude 
continues to increase; however, never reaches the amplitude from the nght side. 

Figure 6.7 shows the FRF between the applied load and the vertical acceleration at 
the pilot heels. The FRF does have symmetry across the tail; however, the magnitude 1s 
substantially less than the lateral. 

Figures 6.8 and 6.9 display the FRF between the same points but at the discrete 
frequencies of 6 and 11.6 Hz. The magnitude was normalized to the maximum amplitude 
FRF. Comparisons can now be more easily drawn for relative values of each FRF. These 
plots indicate the significant role the right side of the horizontal tail plays in transfernng 


energy at both 6 and 11.6 Hz. 
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Figure 6.6: FRF Horizontal Tail to Pilot Heel — Lateral Response. 
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Horizontal Tail - Pilot Heel, Vertical Response 
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Figure 6.7: FRF Horizontal Tail to Pilot Heel — Vertical Response. 
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Figure 6.8: FRF Horizontal Tail to Pilot Heel — Lateral. 
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Figure 6.9: FRF Horizontal Tail to Pilot Heel — Vertical. 
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2. Vertical Tail 


The loading condition and NASTRAN nodes along the quarter chord of the 
vertical tail are shown in Figure 6.10. Again, only one loading condition, loads applied in 
the y direction, was used. As with the horizontal tail, a chordwise investigation was 
conducted. Similarly, the FRF chordwise at five span locations varied only slightly. 


Those results are not presented. 
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Figure 6.10: Vertical Tail Nodes and Loading Condition. 
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Figures 6.11 and 6.12 show the FRF between the vertical stabilizer and the pilot 
heels lateral and vertical. As with the horizontal, no significant response level was evident 
in the 20-25 Hz range. The span axis represents every node along the quarter chord, with 
the base of the tail defined as 0% span and the top representing 100%. 

In Figure 6.11, the lateral FRF of greatest magnitude is at 6 Hz. This occurs at the 
top of the vertical tail. All FRFs decreased as the load was moved down the tail, except 
at 11.6 Hz. The FRF for this frequency reached a maximum at the base of the vertical tail. 

Figure 6.12 shows the FRF between the applied load and vertical acceleration at 
the pilot heels. It is consistent with the information from the previous plot. Magnitudes at 
all but 11.6 Hz decrease as the load was moved down the tail. The amplitudes are 
substantially less than lateral. 

Figures 6.13 and 6.14 display the FRF between the same points but at the discrete 
frequencies of 6 and 11.6 Hz. As with the horizontal stabilizer, the magnitude was 
normalized to the FRF of maximum amplitude, which occurred on the night tip of the 
horizontal stabilizer. 

The top of the vertical has the same magnitude at 6 Hz as does the nght tip of the 
horizontal. The base of the vertical tail has the dominant FRF at 11.6 Hz. Comparing this 
with the pressure transducer data, the nght horizontal tail had the most energy 
concentrated in the 5-8 Hz bandwidth. This is where the 6 Hz FRF is a maximum. For 
the base of the vertical tail, the pressure energy was broadband from 11-19 Hz. The 11.6 


Hz FRF is a maximum at this location. 
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Figure 6.11: FRF Vertical Tail to Pilot Heel — Lateral Response. 
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Figure 6.12: FRF Vertical Tail to Pilot Heel — Vertical Response. 
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Figure 6.13: FRF Vertical Tail to Pilot Heel — Lateral. 
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Figure 6.14: FRF Vertical Tail to Pilot Heel —Vertical. 
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Vil. RESULTS AND ANALYSIS 


A. SOLUTION METHOD 


Having obtained the PSDs of the response and the structural frequency response 
functions, a solution can be obtained for the spectral content of the applied loads. The 
method for accomplishing this is provided by Bendat and Piersol [Ref. 8:p. 240]. A 
summary of this method follows. 

X is defined as a column vector representing the Fourier transform of g input 


records, and Y the Fourier transform of & output records; 


x, y, 

oe ie 
X= Y=|: 

x, LY, | 


X*, Y* = complex conjugate of the vectors X, Y 
X’, Y’ = transpose of the vectors X, Y 


The input and output spectral density matrices can then be defined as 
Zz oe ae : 
G_= in {x x \ input spectral density matrix 


G, = aE fy *y'}= output spectral density matrix 


yy 
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The notation used by Bendat and Piersol for G,, and G,, simply means that the 
PSDs and CSDs of each input and output record are contained with these matrices. 


Another method for representing these terms, and useful for programming is 


Gi, Ge ioe Ge 
o_-|% = es 
Go G. wan Ge 
Ce Cie _ Gry, 
Gy = Oy. Go 5 Oy, 
Siro Gy tee Circe 


Both G,, and Gy are Hermitian matrices, with G,, = G,, for all i andj. Each 
element in both matrices represents either the PSD or CSD for the inputs and responses. 
When the subscripts are the same, as with Gj, the value is a real number and represents 
the PSD for the load at point one. The other elements of the matrix are complex, such as 
G2 which is the CSD between the loads at point one and two. The same is true for the 
response matrix, with the diagonal being the PSDs for responses and the off diagonal 
terms are the CSDs between responses. 


Finally, the FRF, H,,, where x defines the input and y is the output is given by 


ly ly, ly; 

H A 2 2y2 A 2yk 
xy : : : 

ie ees Lage 


Depending on the number of inputs and outputs used in the solution set, the FRF matnx 
may not be square. Using this definition of H,,, it follows that the output Y 1s related to 


the input X by the relationship 


= , fos ° 
Y =H,,X whereH,,, is the transpose matnx to H,, 


Having defined these terms, the following derivation 1s required to arrive at a 


solution for G,,. 


Y = HX 

Nate) 

Yea ia) — tae 
YY =(H/, X")(X'H,,) 


Taking the expected value of both sides and multiplying by (2/T) gives 


G,, =," ¢..H,) 


It follows that 
Gc, =(a',)'¢,(,) 


x xy 


Assuming the required inverse matrices exist, a solution can be obtained for the 
spectral content of the applied loads. In the case of unequal inputs and responses, where 
H,, is not square, the solution, G,x, is found using gaussian elimination. Conditioning of 


the matrices becomes critical during the computations. 


on 


B. MATLAB® PROGRAM 


Using this formulation, a MATLAB” code was written, contained in Appendix E, 
that solved for G,x. First, the user specified the frequency range of interest. This was 
limited from 5-25 Hz because of the range of frequency sweeps use to determine the 
frequency response functions. Then, NASTRAN grid points were selected for application 
of forces. Finally, response points to be used in the solution set were chosen. 

The program looks up, on a frequency by frequency basis, the stored FRFs 
between force application points and response locations. This created the H,, matrix. 
The Gy, matrix was found, again on a frequency by frequency basis, from the selected 
response points. PSDs and CSDs for each set of accelerometer flight test data was stored 
in .mat files. This allowed retrieval for use in the appropnate G,, matrix. Because the 
frequency interval on FRFs was .05 Hz and on the 1024 point FFT PSDs was .293 Hz, 
cubic spline interpolation was used to match the data sets. 

The solution, G,x, represented the PSDs and CSDs of the applied loads at discrete 
frequencies. The frequency band was dependent on what the user had selected for the 
upper and lower limit. The interval between frequencies, or effective bandwidth was 


based on Gy,. Hence for this output, with a 1024 point FFT, the interval was .293 Hz. 
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C. RESULTS 


Numerous cases, with different locations of forces and reponses, were computed. 
Some generalities included all forces placed along the quarter chord of either the 
horizontal or vertical tail. These points are shown in Figures 6.5 and 6.10 and listed in 
Table 7.1. Additionally, the only response data used was from a 100 knot record, entitled 
run 24. The response points contained in this data set are also displayed in Table 7.1. 
Horizontal tail span is defined from the right tip. Vertical tail span begins at the base of 


the tail. 
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Force Application Points 


NASTRAN 
Run 24: Accelerometers Aircraft Location % Span Grid 


Nose Cone Lateral 
Nose Cone Vertical 
Pilot Heel Lateral Horizontal Tail 
Pilot Floor Lateral 


Pilot Floor Vertical 70.70 43008 
Pilot Overhead Lateral : 86.70 4200 
Copilot Overhead Lateral | 80.66 31844 
Vertical Pylon Top Lateral | 65.05 31842 
Hor Stab Left Tip Long 49.45 31840 
Hor Stab Right Tip Long 33.85 31838 
Hor Stab Left Tip Vertical 18.25 31836 
Hor Stab Right Tip Vertical 6.81 30722 
Table 7.1: Response Data and Force Application Points. 
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The first solution set used two forces, one on the night horizontal stabilizer at mid 
semi-span and the other at mid-span of the vertical tail. Then, additional unknown forces 


were included, with different response data. For several frequency bands, the condition 


2° 


number of the FRF matrix, H,,, became exceedingly high. This occurred when the FRFs 
were very small or FRFs from different points were of equal amplitude. In these cases, the 
precision of MATLAB® could not distinguish between numbers and interpreted the two 
rows of the matrix as duplicates of each other. Therefore, the inverse was not reliable. 

Using between two and twelve forces, the response locations were varied to 
produced different solutions. Table 7.2 through 7.8 display the results for seven different 
loading conditions. In all these cases, every response location that run 24 contained was 
used. The tables present the RMS for each force within the frequency bands previously 
identified. In every case, the energy within these three bins represented over 80% of the 
total energy, throughout the 5-25 Hz range, of the force at that node. The final column 
indicates the RMS of the total force. 

The RMS values were determined for each frequency band and then divided by 
the largest RMS value of total force. They represent the percentage of the total force that 
is located at that point on the aircraft and within the frequency range. Similarly, the RMS 
of total force was normalized to this same value. Hence, the RMS of total force that 
equals 1.0 represents the location on the tail that had the largest force throughout the 5-25 
Hz band. Normalizing the forces was done to aid in comparison of loads between cases. 
Additionally, the conditioning of the matrices produced highly inaccurate results within 
certain frequency bands. In these cases, the solution had large amplitude PSDs in 
frequency ranges containing little response and small FRFs. The loads in nine and twelve 
force cases, Table 7.7 and 7.8, are highlighted to show the frequency content and 


magnitudes of the computed loads. 
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Forcing Function - Normalized (lb/Ib 


PE ee ee A Ly ese 


Frequency Bin RMS of 


ae 
| Node | 5-8hz | 10-13hz | 22-25 hz | Total Force 
Table 7.2: RMS of Two Force Solution Set. 















Forcina Function - Normalized (|b/Ib 


je SS ene nn aa en 


rece RM 

| Node | 5-8hz_ | 10-13hz | 22-25 hz | Total Force 

[42008 [| 036 | 042 | 011 | 060 | 
Table 7.3: RMS of Three Force Solution Set. 





Forcing Function - Normalized (Ib/Ib 





Table 7.4: RMS of Four Force Solution Set. 
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Forcing Function - Normalized (lb/Ib 
|| Frequency Bins | «=~ RMS of 
| Node | 5-8hz | 10-13hz | 22-25 hz | Total Force 
I 45108... | 01535 1/0 2S NS | Oo el 0 63 
| 1318449). 0.50. Ono O10 li 060 a) 

Table 7.5: RMS of Five Force Solution Set. 










Forcing Function - Normalized (Ib/Ib 


|r| Frequency Bins ss |:~- RMS of 

| Node | 5-8hz_ | 10-13hz | 22-25 hz | Total Force 

| 42008 | 320.03] 0.01 | 00 | 0 

| 46008 [| 003 | 001 | 000 | 003 | 

| 31844, | 048 | | S| T 

oO a oe iS are 
Table 76. RMS OF Six Force TAG Set. 





RMS of 


. 
| Node | 5-8hz_| 10-13 hz | 22-25 hz | Total Force 


Table 7.7. RMS of Nine Force Solution ser 





1 






Forcing Function - Normalized (|b/Ib 





PCT Frrequencybin |= RMS of 

| Node | 5-8hz_| 10-13 hz | 22-25 hz | Total Force 
0.03813 
| 31842 0.47014 | 0.00037 | 0.53079 
| 31840 | 043749 | 086229 | 0.00037 | 1.00000 






31838 | 0.40227 - 1-0. i re 0.00030 0.86572 

31836 0.17668 |. -0.29943 { 0.00009 0.36249 

30722 0.03467 0.07152 0.00010 0.08449 
Table 7.8: RMS of Twelve Force Solution Set. 








Figure 7.1 through 7.4 show the PSDs of loads for the twelve force solution set. 
The loads on the horizontal tail are within frequency bins which have little response. The 
magnitudes are much smaller than forces on the vertical and associated with poorly 
conditioned H,, matrices. Under the PSD plots is the phase angle between point 31840, 
where the force of largest magnitude is located, and selected points. They were 
determined from CSDs which capture the phasing information. Bendat and Piersol [Ref. 
8:p. 125] provide a method of determining the phase angle from cross-spectral density 
functions. Although the plots are difficult to interpret, if viewed at just the modal 
frequencies earlier identified, 6.0 and 11.6 Hz, the phase relationship can be seen. This is 
available on the vertical tail where forces at these frequencies are located. The horizontal 


tail has no forces at these frequencies. 
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Figure 7.1: PSD of Forces on the Right Honzontal Tail. 
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Normalized PSD - Ib2/hz 
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Figure 7.2: PSD of Forces on the Left Horizontal Tail. 
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PSD of Forces - Vertical Stab (50%-100% Span) 
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Figure 7.3: PSD of Forces on the Vertical Tail — 50% to 100% Span. 
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Figure 7.4: PSD of Forces on the Vertical Tail — 0% to 50% Span. 
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D. ANALYSIS 


Beginning with the initial load set and continuing through twelve forces, the 
computations place the loads of largest magnitude on the vertical tail. Initially, the 10-13 
Hz frequency bin had the greatest energy. As more and more forces were added, the 
magnitude of forces in the 5-8 Hz range grew tremendously. This coincided with a poorly 
conditioned FRF (H,,) matrix at 5 Hz. Even when eliminating the contribution to RMS 
from this frequency, forces from the base of the vertical tail to mid-span, at both 5-8 and 
10-13 Hz, play the dominant role in producing the response. 

In going from nine to twelve forces, the 10-13 Hz frequency bin, again, became the 
dominant source of energy. The problems with the conditioning of the FRF matrix, 
discussed above, were again present in the final solutions. The point where the largest 
forces were placed was grid 31840. Figures 7.5 and 7.6 display the distribution of forces 
along the horizontal and vertical tail. It can be seen that more energy was placed on the 
lower two-thirds of the vertical tail in both the 5-8 and 10-13 Hz frequency band. This is 
consistent with the results from the nine load case, shown in Table 7.7. 

The results also match the information from the analysis of the frequency response 
functions and pressure transducers. The solution places loads on the fuselage at points 
that have high magnitude FRFs at 11.6 Hz. Pressure transducers confirm broadband 
excitation, including this frequency, and of largest amplitude at the base of the vertical tail. 

Although the solution indicates a force on the vertical tail in the 5-8 Hz frequency 
bin, pressure transducers do not support this. The solution also has no forces in the 22-25 


Hz range. The lack of NASTRAN modes in this bin has previously been discussed. 
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Figure 7.5: Distribution of Forces, 12 Force Solution Set - Horizontal Tail. 
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Figure 7.6: Distribution of Forces, 12 Force Solution Set — Vertical Tail. 
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Vill. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


The goal of this work was to quantify the magnitude and location of aerodynamic 
loads contributing to the random excitation of the empennage. Then, follow-on efforts 
would have the proper focus and direction in modifying the structure. This was 
accomplished by first conducting an extensive analysis of flight test response data and 
pressure transducers. Due to the nature of the accelerometer time histories, spectral 
analysis was used. With this accomplished, the mode shapes and frequency response 
functions of the aircraft were determined. A NASTRAN model, modified to replicate the 
test flight aircraft, provided this information. Finally, a MATLAB?® program was written 
to compute the spectral content of the aerodynamic loads producing the random 
excitation. 

Each portion of the analysis provided insight into the nature of the response, 
locations of excitations, and the method in which the structure transfers input to output. 
From all of these elements, combined, conclusions were drawn. 

The dominant response is lateral and occurs in three frequency bands or bins. The 
first bin is 5-8 Hz, then 10-13 Hz, and finally 22-25 Hz. Each accelerometer has its 
maximum value within one of these bins. The top of the vertical pylon has maximum 
response at 6.00 Hz. For the nose cone, pilot’s station, and horizontal tail, 10-13 Hz 
excitation is the largest. Finally, the 22-25 Hz response is restricted primarily to the CPG 


station and vertical acceleration of the horizontal tail. 


11] 


Within the first two bins, there are mode shapes that correspond with the motion 
seen in the output of the 22 accelerometers. The aft fuselage lateral torsion mode at 6.00 
Hz and horizontal stabilizer yaw mode at 11.6 were identified not only through 
NASTRAN but also during actual frequency response hover testing. No mode was 
identified or significant frequency response function found within the 22-25 Hz bin. The 
horizontal stabilizer vertical bending mode at 26.02 Hz matches the response seen in the 
22-25 Hz bin. The actual natural frequency of this mode was not identified during flight 
test; however, in previous NATRAN models it is at 22.06 Hz. 

From the forces computed through the MATLAB” program, good correlation 
existed with the previous analysis. Numerous loading combinations were applied with 
maximum loads continuously placed on the vertical stabilizer from the base to 50% span. 
The 10-13 Hz frequency bin contained the largest force. This agreed with data from 
differential pressure transducers which showed the highest broadband pressures at the base 
of the vertical tail. 

From this, it can be concluded that forces within the 10-13 Hz bandwidth 
impinging on the lower two-thirds of the vertical tail play the most significant role in 
creating the random excitation. This is supported by the differential pressure transducers. 
Solutions want to place 5-8 Hz forces in this same area. This is not supported by pressure 
transducers and leads to the conclusion that main rotor 1P 1s the primary source of 6 Hz 
excitation. Because the 22-25 Hz excitation is so localized and no FRFs of magnitude are 
in this region, the least squares solution provides no contributing loads. If the horizontal 


stabilizer vertical bending mode is not at 26.02 Hz but lower, the broadband excitation 


Zz 


seen in the differential pressure transducers can excite this mode and produce the recorded 


response. 


B. RECOMMENDATIONS 
1. Correlate with Wind Tunnel Testing 


Boeing and Sikorsky have conducted two wind tunnel tests. The purpose was to 
quantify the flow impinging on the tail and identify changes to the surface of the aircraft, 
such as strakes, that would modify the flow and reduce excitation levels. The results of 
this analysis should be compared to data from the wind tunnel. This could confirm the 
conclusions of this study about the nature of the aerodynamic loads and their primary 
locations on the fuselage. 


2. Parameter Identification and Correlation with NASTRAN 


During hover frequency response testing, actual natural frequencies were identified 
for three modes. As additional test flight hours become available, more parameter 
identification can occur. This will better identify the exact frequency of the horizontal 
stabilizer vertical bending mode, as well as clarify other modes. In addition, the change in 
natural frequencies due to adding stiffening brackets and tip weights to the horizontal 
stabilizer can be assessed. With this information, a detailed analysis of the NASTRAN 
model can be conducted with modification of natural frequencies and mode shapes to 


replicate the actual aircraft. 
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3: Additional Airspeed Analysis 


The computed solution using the method of Bendat and Piersol was done for only 
100 knot records. After receiving additional airspeed records from Sikorsky, the PSDs 
and CSDs were computed; however, no complete solutions were obtained for the spectral 
content of forces. Calculating solutions throughout the range of airspeeds can provide 
further confirm and clarity as to when the random excitation begins and if the 6 Hz 
response is due to main rotor 1P. 


4. Structural Optimization 


Ultimately, the goal is to eliminate or reduce the response to acceptable levels. 
Using the conclusions from this analysis and confirmation from the sources discussed 
above, the structure can be altered to reduce the frequency response functions within the 
excitation bandwidth. This can be accomplished using several different methods, or, 
optimally, a combination of all. 

First, using the NASTRAN model, structural modifications can be done that move 
the natural frequency, 11.6 Hz, of the horizontal stabilizer yaw mode. The focus of the 
changes should occur on the lower two-thirds of the vertical tail. Extensive work was 
done by Tobin [Ref. 12] and Shoop [Ref. 13] to increase torsional stiffness of the aft 
tailcone. The work attempted to stiffen the tail and move the aft fuselage lateral torsional 
mode away from MR 1P. Their analysis and strain energy diagrams provide insight to 
methods of altering the structure, specifically for the 11.6 Hz mode. 

Second, changes are proposed to the design of the tail fold, spars of the vertical 


tail, and the tail fan gear box struts. Additionally, a trade study is being conducted by the 


114 


Comanche Program Management Office (PMO — Comanche) on the mounting of the tail 
landing gear bay. As these designs are finalized, the dynamics of the tail will change. The 
quantification of the flow determined through this analysis should aid in creating a 
dynamically “tuned” design. 

Finally, modifying the flow around the main rotor hub and pylon has been 
investigated through the wind tunnel experiments previously discussed. Modifications to 
the structure must be done in conjunction with changes to the flow. Only through this 
integrated approach can an acceptable solution be obtained. In this manner, both the flow 
and structural response can be optimized for maximum reduction of the airframe vibratory 


response. 
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APPENDIX A: TEST FLIGHT PARAMETERS 








| Affected Components | Units _ 


Horizontal Stabilizer find 


\Parameter 


ae ne Gbneer 
A ae Bending, LHS (near Root 


RMR Horz Stabilizer Flatwise 
aE Bending, RHS (near Root 





Horizontal Stabilizer 






Horz Stabilizer Flatwise - 

HBNORML2 Ben nLIst onitoard Horizontal Stabilizer 
Horz Stabilizer Flatwise . Bs 

HBNORMR2 Econ erc\(cuthoard Horizontal Stabilizer 





HSTBFWDpB | Hor Stab Fwd Folding Joint Horizontal Stabilizer 
Bending LHS 


HSVAFTL | Hore Stab ARLHS Spar Shear | 544+ ontal Stabilizer 
Force — lb 
HSVAFTR aang Aft RHS Spar Shear | porizontal Stabilizer 


HSVFWDL_ | Ho! Stab Forward LHS Spar | by) ontal Stabilizer | N/A 
Shear Force — Ib 


HSVEWpR _ | Ho! Stab Forward RHS Spar | proc ontal Stabilizer | N/A 
Shear Force — | 


b 

VSFFTS1 Vert Stab Attach Fwd Ftg Outer Vere eion ae? 
Cap LHS 

VSFFTS?2 Vert Stab Attach Fwd Ftg Outer A ereee sone “wy 
Cap RHS 

VSAFTSI Vert Stab Attach Aft Ftg Outer Vert / Shroud a“? 
Cap LHS 

VSAFTS?2 Vert Stab Attach Aft Ftg Outer Tchad au“? 
Cap RHS 


Table A.1: Test Flight Parameters. 







j= 








es 












‘Mnemonic | Parameter — 8 | Affected Components — 


BVSRTMX | Vertical Stab RG ot Roll Rehan d oe 
Moment 

DHSDRAG Horizontal Stab Drag - DP Horizontal Stab | N/A | 

DHSTABMX | Net Stab Roll Moment -DP _| Horizontal Stab ) inlb 


DHYAWMZ Honizontal to Vertical Stab ee cea 
Yaw Moment — DP 


DHSTABVZ_| Net Stabilizer Lift (Ib)- DP _| Horizontal Stab 


DVSAFTLH Derived Vert Load Aft LHS Sever 
Lug - DP 
DVSAFTRH Derived Vert Load Aft RHS er hennt 
Lug - DP 
DVSFWDLH a ceases Load ae ment cinoud eal 
DVSFWDRH mop’ ta LATE) SS | op a) Sto a 
aie Stab Root Roll MMT | Vertical Stab 
a Stab Root Pitch Vertical Stab 
DVSTABVY Derived Vertical Stab Lateral Vert / Shroud 
Force — DP 
HBEDGELI Horz Stab Edgewise Bending, oon sen 
LHS (near root 
HBEDGERT © | ueeeee eee Seen ie rior onnsTat 
RHS (near root 


Table A.2: Test Flight Parameters. 
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APPENDIX B: PSD AND CSD PROGRAMS 


% psaetltrdata.m 


% Program reads the input flight data for a run and the computes the 
% PSD for each accelerometer time history 

% The input record is an 18,000 X n matrix with each column 

% representing another transducer 

3 Input the data set you want to read, sampling rate, data set you 


3 want plotted 
clear 


Pune Labbe Ws (tile 


dito 417 300; 

string =['runl40']; 

speed =['140 Knots']; 

eval (['lload™', string, ';"1); 
i = 0; 


for col = 4332 
= ee 


n = 1024; 
noverlap = round(.67*n); 
eval({'x = *,string,' data({:,col);‘])); 


$ Using the SPEC command 
$ [Pxx,F] = spec(x,n,noverlap,hanning(n),1/dt); 


$ Using the PSD command 
[Prec ep sa(x,n, 1/dt,hanning (n),noverlap) ; 


$ Normalize by multiplying by 2/frequency sampling rate 
$ To get the rms multiply by sampling rate/nfft then add 
$ This is from Pareseval's theorem 


Pee (1 / de) jee PxXx; 
% Another way to normalize so that the rms is just amplitudes add 
$ together 


$ With this method, the value from different length fft's is the same 
SK — ae el, ite) aS O0/n)  * Pxx; 


119 


Soe OLEInG Punecrion 


plot _ max = round(25/((1/dt)/n)); 
ymax = max(Pxx(1l0:plot max,1)); 
figure (i) 
Plec(F/Pxx(s, 1)) > 
@pieatel - 
axis([0 25 0 ymax]); 
set (Gea, “xXtick 7 10rieZoi 
eval ([‘title(*' PSD ', deblank (eine rabets (comm) 75 +) DUMZSEr (aye borne 
Fors) ee, 
xlabel ("Frequency - Hz'); 
ylabel('g>s*2/Hz'); 
eval ( [* text(5,.9* ymax,'* Run *, speed; 
', Hanning Window, 67% Overlap'')']); 


% Save the PSD data in another matrix for use in solving for forcing 
fLunceLons 


* ‘Name Gf the matrix eseved 125 psa runza eimaee Onmar mand has che 
frequency and 


% the psd for all accelerometer data we had for run24. 
eval | psa. ,string, a3, 2) = Bsoe(s, Pye, 
end 


$save frequency as a separate variable 
pSdmiabets = run iabels(4-size( rune vabela hy 


eval(@i save psd “,string,” psde, String, f psce rebels) ji, 
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% GsGerie Mata. m 

% Program reads the input flight data for a run and the computes the 
% CSD (cross spectral density) between output accelerometer time 

% hiskeory 

% The input record is an 18,000 X n matrix with each column 

% representing another transducer 

% Input the data set you want to read, sampling rate, data set you 


% want plotted 


at = 17300; 

load run24 

k = 0; 

fon, col = soc. 

for 2 = Lsi2 

k =k4#+i1; 

fobs oS 

n = 1024; 

noverlap = round(.67*n); 


x = run24 data(:,col); 
y run24 data(:,1i1i); 


% Using the CSD command 
fPxy, Sle=eesdi(x, y,a,1/at, hanning (n) ,noverlap); 


% Normalize by multiplying by 2/frequency sampling rate 
% To get the rms multiply by sampling rate/nfft then add 
%$ This is from Pareseval's theorem 


Pxy = (2/(1/dt)) * Pxy; 


%$ Another way to normalize so that the rms is just amplitudes added 


% together 
% With this method, the value from different length fft's is the same 


weve ey lediz)) = (3007/n) * Pxy; 


SePeLOcEINg FuNncr lon 

plot _min = round(5/((1/dt)/n)); 

plot _ max = round(25/((1/dt)/n)); 

ymax = max (abs (Pxy(plot_min:plot_max,1))); 
figure (k) 

Slo Cie a,b, dna, 

Bloe (habs (Pxyvii, i) 

ogslic |= 

axis([5 25 0 ymax])); 

Set(gcea, xtiek', (9:17:25) ).- 
Svaliecrtlcel espe: ,deblankinun24 labels (col,:)), = =, 
deblank(run24 labels(ii,:)),' °')°"]); 
$xlabel('Frequency -— Hz"); 
ylabel('g’s*2/Hz (Amplitude)'); 
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Subpmoe (2,2) 
plot (F,180/pie* anglie( Pxyieel) )); 
angle max = max(180/pi * angle(Pxy(plot_min:plot_max,1))); 
arteiion2co lee Leu |); 
Seeigea, “Xtick | lo+l:25)); 
set (gea, ‘ytick', |-l8O0s26G2isi0)))) ; 
efigal.e | 
xXlabely "Prequency .— Jhz. ), 
ylabel('Phase Angle - degrees') 
eval ([’ title” "Run “7numZser{ runZagedaca( il). 
";, Hanning Window, 6/2 "Overlap, . ~;numzstr(n), Point FER), 
seval ({'text(S,.89*angle max, '*Phase Angle 
',numéstr(run24 data(1,1)),... 


% ', Hanning Window, Gis Overlap, @mum-str(n)>.~ Point. FFE saa, 
$ Save the CSD data in another matrix for use in solving for forcing 
%$ functions 

%$ Name of the matrix saved is psd_run24 (mat format) and has the 

$ frequency and the psd for all accelerometer data we had for run24. 
eVab | escecune4  pmumZser (COM se) ee ep 7 ae, 

end 

seval({| "save psd run24 esd runz4 7 numzscer(col—3), fF —appenda, 1); 
end 
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APPENDIX C: MODE SHAPES 


Mode shapes for selected modes are presented on the next eleven pages. 
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APPENDIX D: ACCELEROMETER LOCATIONS & NASTRAN GRIDS 


RAH66- MODEL J yrahj194.dat 4/17/98 sp 
COMPARISON OF KEY FLIGHT TEST LOCATIONS WITH SUGGESTED NASTRAN APPROXIMATE LOCATIONS 


FLIGHT TEST NASTRAN 

Approximate Location Location 

ITEM [LOCATION STA BL WL STA BL 

No mm mm mm mm mm 

VNOSE:D:300 { 5286.00 000 2720.00} 5298.00 0.00 
LNOSE:D:300 
LONOSE:D:300 


PILOT HEEL LPHEEL:D300 
VPHEEL:D300 


VCOPHEEL:D:300 


VCOPF:D:300 


LPOVRHD:D:300 | 7420.00 


COPILOT OVERHEAD LCPOVRHD:D:300] 8795.00 


VERTICAL PYLON TOP LVPYLNT:D300 


HOR. STABILIZER LEFT TIP VSTABTLT: D:300 
LOSTABTL:D:300 


HOR. STABILIZER RIGHT TIP VSTABTRT:D:300 
LOSTABTR:D:300 


HOR. STABILIZER VERTICAL HBNORML1:D:300 

BEND. LT SIDE 
-281.41 
-281.41 





Table D.1: Accelerometer - NASTRAN Locations. 
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FLIGHT TEST 
Approximate Location 
ITEM JLOCATION BL WL 
No. 


HOR. STABILIZER VERTICAL HBNORML2:D:300 
BEND. LT SIDE 


HOR. STABILIZER VERTICAL HBNORMR1:D:300 
BEND. RT SIDE 


HOR. STABILIZER VERTICAL HBNORMR2:D:300 
BEND. LT SIDE 


HOR. STABILIZER EDGE HBEDGEL1:D:300 
BEND. LT SIDE 


HOR. STABILIZER EDGE HBEDGER1:D:300 
BEND. RT SIDE 


LT AFT SHROUDIV. TAIL VSAFTS1:D:300 
INTERFACE FTTING 


RT AFT SHROUD/IV. TAIL VSAFTS2:D:300 
INTERFACE FITTING 


LT FWD SHROUDIV. TAIL VSFFTS1:D:300 
INTERFACE FTTING 


LT AFT SHROUDWN. TAIL VSFFTS2:D:300 
INTERFACE FTTING 


FWD H. STAB/V.TAIL - DHSTABMX: D300 
FITTING CL DHSTABMZ:D300 
DHSTABVZ:D300 





Table D.2: Accelerometer - NASTRAN Locations. 
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APPENDIX E: LOAD CALCULATIONS PROGRAM 


% ff solve nonsym inter.m 

% Program that solves for the forcing functions, in PSD by using the 
% techniques discussed by Bendat and Pierson for multiple input 

% single output 

% Gx oi |) le = Gyy * 1 (Hew) 7-2 

% First must identify the location (nodes) of the forcing function 

% the location of the response (nodes), and the direction of the 

% response 

% The directions are somewhat backwards - due to the way I had Nastran 
$ write out the punch files ( which became .mat files) 

$ x = 3 (longitudinal) 

$ y = 2 (lateral) 

$ ze=t1 (vertical) 

% The response data uses the following script: 


$ First is node in Nastran Mode.1l, next is direction and last is the 


$ pointer into the psd matrix (which column) 
nel = [ele 41 2 1}; % Nose Cone Lateral 
phl = fLGa0s 2 2}; % Pilot Heel Lateral 
pfl = [7246 2 3}; %& Pilot Floor Lateral 
pol = [7514 Z 4}; % Pilot Overhead Lat 
cpol = [8614 2 5}; % CPG Overhead Lat 
ncv = [o1e41 1 6}; % Nose Cone Vertical 
piv = [7246 1 jis Sele Elocr Vertrcal 
vptl = [32341 2 8); % Vert Pylon Top Lat 
hrtv = [48108 1 913" 3 HerStab Rt Tip Vert 
hitv = [41108 1 LO} =e HOrotab be lip VenE 
hrtlo = [48108 3 Ll ie ceOrcS cab Re Tip long 
hltlo = [41108 3 iZ\ 7" “seenoroctao Lt fap. Long 
% Use the following format to enter the input-ouput locations and 


% 


run 


specify which run set 


= 24; 


sno fit = 4096; 


ee 
$ x 
% 

fore 


orce = (326 46008 } 
esponse = [ pol 
phl]; 
e = [ 42008 44008 46008 31844 31840 31836 ]; 
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response = [ phl 


piv 
Cool: 
mer 
ncv 
Vie 
hie Vv. 
heey 
betlo 
hlLeLoo 
Assumptions is that loadsi are correlation - are the loads correlated? 
If the loads are completely uncorrelated then set corr = 0 
If it is something inbetween then define the matrix relate 
force Cotk a= Ay, 
RESPONSE COmr cs ain, 
Use the following format for relate - forcing function(load) to load 
(which are correlated) 
tones relate — ) (bev 
> i 
001) This tells us force 1 is related to force 2, 
all others are not 
force melate = [iv 
Ol); 
Load the frequency response function matrix 
Counter to use for matrix index 
= 0; 
[Neer Ow, tOrce num) = size(force); 
[tesponse num,ne col] = size(response); 


BOG f eessce Vy) 
$ frf file must have the name format frfnode 


For) © — i:response num 
k = k#+1; 
eval —' toad f£rl' snumZser(torce (1, f) )* se ai 
CV at ence ies Ere 7 numZscte (response (i i) aI Ri 
Vow (iret sel (jaa reer setts ae (fpEssekicw) + 
PEE SSS 1 C.74 )aeey) (Let Sel seed 2 et Seinen ee D8) lu, 
H(:,k) = y(:,response(r,2)); 
end , 


end 
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H (frf matrix has now been loaded in the following format; from our 


g 
$ example nodes - 

See = Prete? Mis HZ h22en2s HSlsH3s2 H33]; 
% 


$ Because of the weird way I saved the pch files, the vector {time} nS 
> Seactlally ~Lrequency 


Now we need to load the response data from the PSD and CSD we have 
% done before 


oe 


eval (['load DsdacUpe, Mune Serd Eun) y, oN) 


eee Se eon Gee) 


eval(['load psd_run',num2str(run),'_',num2str(no_ fft),'; ‘']) 
end 
a 0; 
for £ = li: response num 
forwr = li responsesnum 
k = k+l; 
if f == 
G(:,k) = psd _run24(:,response(r,3)); 
else 
eval(['G(:,k) = ... 
esd _run24 ',num2str(response(f,3)),'(:,response(r,3)); baa) 
end 
end 


end 


% Form the Hxy and Gyy matrix for each individual frequency, 


® solve for Gxx 
% and save in another matrix 
kk = QO; 


Sere ree oy 
end freq = 24.00; 


PEmCnE teen em ySiea ine |= min(abs (start. freq - F)); 
[i ditt; ieindex end] = min(abs (end freq — F)ge 


fom tte -  feindex start... 9 ei ndexgend 
Geccunes— OF 
k = QO; 
kk = kk + 1; 
Pei wee sponse sui 


for Miss). sesponse num 
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cP CP ® &H 


oP o® o® oP 


Load the Gyy square matrix 

Gyy (7m) GEE, Geeoune snipe 

end 

G count = G count + response_num,; 
end 


For 
for m 


f toree num 
l:response_num 


Load the Hxy square matrix 

[Mindi Evel sandex | —*mbmlacs (a (rr > cima aie 
Cubic spline Ineerpola een 

Hxy (2m) = interpol (time, H(i, ktm), bitty spline’); 
Linear Interpolation 

Hxy(l,m) = interpl(time,H (3, kim) ee (tt i spline 
end 


kK = K t+ Fesponse num; 
end 


tf extscmemeerte late ja=— 1 
Hxy = force relate. =~ Hxy, 
end 


eee xe Geo Conse mre ictus aie 
Gyy = excite wrelate ~*~ Gyy, 
end 


ape ae ceysleoke = Mele icig 
Gyy = eye(size(Gyy)) .* Gyy; 
end 


Now we have to solve for the Gxx matrix via the technique 
discussed above 


ie mormee weOra) —— 0 NcetrcoDOn oc mCere=e a — at) 
Gyy = diag(Gyy) ; 

Gxx = ((abs(Hxy')).%*2) \ Gyy; 

else 

Gxx = Hxy'\ Gyy /Hxy; 

end 


Saving the output data Gxx to a matrix of the following form 
Gxx = [Gll Gl2 G13 G21 G22 G23 G31 G32 G33) 

Also outputs a vector, Freq, that has the frequencies 

each row represents 


(Gxx row, Gxxecol) "= size(Gx=)7, 


if Gxx col == 
G output(kk,:) = Gxx'; 


else 
k = 0; 
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EGG) al Gx EOw 
for mm = 1:Gxx_col 
G output (kk,k+mm) = Gxx(1ll,mm); 
end 
k = k+ force num; 
end 
end 
frequency(kk) = F(ff); 
cond Hxy(kk) = cond(Hxy); 
rank _Hxy(kk) = rank (Hxy) ; 
end 
big “Kalantar Pesequeney ~GrourpuLr cond Hxy' rank _Hxy']; 
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